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The MiniBooNE experiment is a — > v e and P M — > P e appearance neutrino os- 
cillation experiment at Fermilab. The neutrino mode oscillation analysis shows 
an excess of u e candidate events in the low-energy region. These events are an- 
alyzed under the SME formalism, utilizing the short baseline approximation. 
The preliminary result shows the time independent solution is favored. The 
relationship with the SME parameters extracted from the LSND experiment 
is discussed. The systematic error analysis and antineutrino mode analysis are 
outlined. 



1. MiniBooNE low-energy excess 

The main goal of the MiniBooNE experiment 1 at Fermilab is to confirm 
or reject the LSND oscillation signal. 2 The LSND experiment claimed a 
v e appearance signal from beam by fi + decay at rest (~40 MeV), and 
this corresponds to Am LSND ~ O.f — f.O eV 2 . The MiniBooNE oscilla- 
tion analysis in neutrino mode rejected this Am 2 SJV£ , 3,a because Mini- 
BooNE did not see the excess in the energy region where the LSND signal 
is expected {E® E > 475 MeV) under two-neutrino massive model. How- 
ever the first oscillation result had unexplained excess at low energy region 
(EQ E < 475 MeV). This was confirmed after a year long reanalysis. 5 This 
signal cannot be understood from the three-neutrino massive model, but can 
be understood with a Lorentz violating neutrino oscillation model. Espe- 
cially, a neutrino oscillation model based on Lorentz violation, the so called 
tandem model, 6 predicted the low energy excess of MiniBooNE. Therefore, 
it is interesting to search for sidereal variation in these low energy excess 
events. 



a A recent analysis shows an excess for the antineutrino mode. 4 
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2. Data set 

The data set used in this analysis is limited to the neutrino mode. Also, all 
the results presented here are preliminary. 

We used neutrino mode data from March 2003 to January 2006, and 
October 2007 to April 2008, corresponding to total of 6.361 x 10 20 protons on 
target (POT). In this data set, we found 544 v e candidate events in the low 
energy region (200 MeV< E® E <475 MeV), and 420 v e candidate events 
in the oscillation candidate energy region (475 MeV< E® E <1300 MeV). 

Before we proceed to the extraction of Standard-Model Extension 
(SME) parameters from our data set, we performed statistical tests to quan- 
tify any statistically significant time variations. We employed Peason's x 2 
(P"X 2 ) test and the Kolmogorov-Smirnov (KS) test. 7 The results are shown 
in Table 1. For the P-x 2 test, we chose the bin size so that the predicted 
number of events in each bin is greater than 5. The KS test is performed on 
the unbinned data set to maximize the statistical power. There are 4 data 
sets, the low energy data sample and the oscillation candidate energy data 
sample, each of which include GM time distribution and sidereal time distri- 
bution parts. None of the 4 data sets show statistically significant variation, 
in other words, all data sets are consistent with the flat hypothesis. 

Table 1. A summary of preliminary results from statistical tests on the 
sidereal and GM time distributions of the v e candidate data. The null 
hypothesis tests compare the data with a constant time distribution. 

null hypothesis tests (preliminary) 
low energy region oscillation candidate energy region 
sidereal GM sidereal GM 



# of events 544 






420 


Pearson's \ 2 '- 








AWs " 107 


107 


83 


82 


X 2 107.6 


106.0 


69.6 


76.2 


P(X 2 ) 0.47 


0.51 


0.85 


0.66 


Kolmogorov-Smirnov: 






P(KS) 0.42 


0.13 


0.81 


0.64 



3. Analysis and results 

Although there is no statistically significant time variation in any of the 
data sets, small time variations are not rejected. So we continue to search 
for Lorentz violation in our data sets. This analysis follows one performed by 
LSND collaboration. 8 We use the SME formalism in the neutrino sector. 9 
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The effective Hamiltonian written under the SME formalism is used to 
write down the neutrino oscillation formula, with the assumption that the 
neutrino baseline is short enough compared to the oscillation length (short 
baseline approximation) . Then, we can get the following oscillation formula 
for the to v e appearance signal, 10 



Here L is the baseline of neutrinos, Cl<g is the sidereal frequency (=2-7r/23h 
56min 4.1s), and T ffi is the sidereal time. The 5 parameters (C) eM , {A s )en> 
{•Ac) eiJ ., {B s ) etl , and (£> c ) eA1 depend on the SME coefficients (az,)" and 



We use the unbinned likelihood method to find these parameters. How- 
ever, in these proceedings, fitting parameters are limited to the first 3 
parameters {{C) ell , {A s ) et i, and {A c ) e n)- This 3 parameter fit can be in- 
terpreted as the case where Nature only has CPT-odd Lorentz violating 
coefficients. 

The fitting results for the low energy region data samples are shown 
in Fig. 1. The solutions are duplicated due to the nature of Eq. (1). From 
Fig. 1, one can see that only {C) ell (sidereal time independent term) dif- 
fers from zero with statistical significance. This means that the data favor 
the flat solution, and this is consistent with the flatness tests in previous 
section. Also, the increase of goodness-of-fit after the unbinned likelihood 
fit is small (P(x 2 ) = 0.60 ->■ P(x 2 ) = 0.77). This is expected because the 
flat hypothesis already provides a very good fit. The situation is slightly 
different for the oscillation candidate energy region. Since the excess itself 
is small, all parameters are statistically consistent with zero. 

Finally, the SME coefficients with la statistical errors are extrapolated 
from the fit results. Note, 350 MeV is used as an averaged neutrino en- 
ergy for the low energy region data sample. From (C) eAi , an extrapolation 
is performed assuming only one SME coefficient is nonzero. {A s )en and 
{A c )efi are solved assuming two coefficients are nonzero. The results are 
summarized in Table 2. We also calculate SME coefficients from the LSND 
experiment. 8 The 3 parameter fit result from LSND has 2 solutions in la 
area, but we accept the statistically significant solution. We used 40 MeV 
for the averaged neutrino energy. 



For the MiniBooNE result, we can see (ol)^, (ol)^, (cz,)fj, (cj,)^f , 
and {cl)cu are statistically significant from zero. This is the consequence 



P 



7^-72 ( C W + {As) e» sincj ffi T e + {Ac) en COSW ffi T e 

+ (B s ) eA1 sin2w e T e + (B c ) efl cos2uj (d T (B | 2 . 




(1) 
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(As)^, lowE MB sideal time (sec) 



Fig. 1. The preliminary unbinned likelihood fit solution for the low energy region side- 
real time distribution. The first 3 plots show parameter space with the preliminary best 
fit point and lcr contours, and the last plot shows the preliminary best fit solution and its 
goodness-of-fit. The curve shows the model line with the preliminary best fit parameters, 
whereas dotted and dashed lines show the flat solution and background. The upper x 2 
and P(x 2 ) are after the fit, the lower ones are before it (flat solution goodness-of-fit). 



of (C) efl being nonzero from the 3 parameter fit result. On the other hand, 
since (*4 s ) g/J is nonzero for LSND, (a L )*^ (cl)s* , and (c L )~f are sta- 
tistically significant. Since the MiniBooNE result provides neutrino SME 
coefficients, and the LSND result provides antineutrino SME coefficients, 
these two results are not inconsistent. 

In conclusion, we performed a search for Lorentz violation in the Mini- 
BooNE low energy excess sample. We discovered that the time independent 
solution is favored. To complete the analysis, we are planning a study of 
systematic errors. We also plan to extend this analysis to the antineutrino 
mode. This is very interesting because the antineutrino mode has an excess 
in the oscillation candidate energy region. 4 
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Table 2. The preliminary extrapolated SME coefficients from the 3 parameter fit 
result. The 3 missing rows corresponds to (S s )e^ an d (Bc)e^, which we did not use 
for the fit at this time. We use 350 MeV (MiniBooNE) and 40 MeV (LSND) for 
the averaged neutrino energies for these extrapolations. The numerical entries are 
preliminary extracted values with lcr statistical errors. 



Preliminary SME coefficients from the 3 parameter unbinncd likelihood fit 



SME 
coefficient 



result from 
MiniBooNE 



SME 
coefficient 



result from 
LSND 



(°-lYL 
<Al)L 

{c L )fj 
(cl£ z 

(0L)f Z 



(-3.1 ±0.9) x 10- 20 GeV 
(0.6 ± 1.9) x 10- 20 GeV 
(-0.9 ± 1.8) x 10- 20 GeV 
(-4.2 ± 1.2) x 10- 20 GeV 
(7.2 ± 2.1) x 10- 20 
(-0.9 ± 2.8) x 10" 
(1.3 ±2.6) x 10" 
(5.9 ± 1.7) x 10" 



-j-20 

-20 
-20 



(-1.1 ±3.7) x 10- 20 



(1.7 ±3.4) x 10 
(2.6 ±0.8) x 10 



-20 

-19 



(ol)|z 
{clY, 



(cl) 
(cl) 
(cl 
(cl) 
(cl) 
(cl) 
(cl)\ 
(cL)fn Z 
(cl)I Z 



eu 
TY 

>P 

eu 
,XY 



,X_Z 



(0.2 ±1.0) x 10- 19 GeV 
(4.2 ±1.5) x 10- 19 GeV 
(-1.7 ± 1.8) x 10~ 19 GeV 
(1.0 ±5.4) x 10- 19 GeV 
(0.3 ± 1.8) x 10- 18 
(-5.2 ± 1.9) x 10~ 
(2.1 ± 2.2) x 10~ 18 
(1.3 ±6.7) x 10- 18 



-18 
18 



(-2.7 ± 1.0) x 10" 

(1.1 ± 1.2) x 10" 
(-1.1 ± 5.9) x 10 



17 



-17 
-18 
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